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1.  Introduction 

Brittle  solids  such  as  ceramics,  rocks,  minerals,  and  glasses 
demonstrate  a  variety  of  deformation  mechanisms  when  subjected 
to  shock  loading  at  stress  levels  above  their  Hugoniot  Elastic  Limit 
(HEL,  i.e.,  yield  strength  under  planar  shock  compression).  In 
crystalline  materials,  these  mechanisms  can  include  intergranular 
fracture,  transgranular  fracture,  pore  collapse,  dislocation  glide, 
stacking  fault  propagation,  twinning,  phase  transformations,  and 
shear  localization.  Brittle  materials  are  informally  distinguished 
from  ductile  solids  such  as  metals  by  their  increased  tendency  to 
fracture  rather  than  deform  plastically  by  slip  or  twinning.  Bonding 
tends  to  be  covalent  or  ionic  in  character,  and  though  exceptions 
exist,  brittle  solids  tend  to  have  a  relatively  large  ratio  of  shear  to 
bulk  modulus  (small  Poisson's  ratio)  in  comparison  to  ductile  solids 
[1]. 

In  the  context  of  dynamic  loading,  engineering  ceramics  tend  to 
have  a  large  HEL  and  small  spall  strength  relative  to  engineering 
metals.  In  brittle  solids,  dynamic  yield  strength  also  tends  to 
depend  more  strongly  on  pressure  [2,3].  Physically,  such  pressure 
dependence  results  from  frictional  resistance  to  crack  sliding  and 
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resistance  to  dilatation  with  pressure  that  accompanies  such 
sliding.  Pressure  dependence  of  strength  varies  among  brittle  solids 
with  different  compositions  and  microstructures  [4,5].  As  shock 
pressure  increases  beyond  the  HEL,  shear  strength  may  increase, 
decrease,  or  remain  constant  depending  on  the  material  [6],  but 
spall  strength  more  often  degrades  as  cracking  increases  with 
increasing  shock  pressure  [4,7],  in  the  limit  the  pulverized  material 
having  zero  tensile  strength.  In  contrast,  in  ductile  solids  deforming 
by  dislocation-based  slip,  the  dependence  of  shear  strength  on 
pressure  is  normally  small  (e.g.,  Peierls  barrier  depending  on 
pressure  through  the  shear  modulus  [8]),  as  is  the  dependence  of 
spall  strength  on  impact  stress,  though  again  exceptions  are 
possible. 

The  first  part  of  this  paper  (§2)  develops  a  finite  strain  model  for 
brittle  solids.  A  multiplicative  decomposition  of  the  deformation 
gradient  is  used,  wherein  inelastic  deformation  arises  primarily 
from  micro-cracks,  though  contributions  from  other  mechanisms 
such  as  dislocation  motion  are  not  excluded.  This  geometrically 
nonlinear  approach  parallels  that  of  finite  elastic-plastic  theory 
typically  used  for  ductile  solids  [8-10]  and  differs  from  rate-type 
approaches  more  often  encountered  for  computational  modeling 
of  brittle  solids  under  dynamic  loading  [11—13].  The  present  work 
focuses  on  an  internal  energy-based  thermodynamic  formulation, 
with  internal  energy  depending  on  entropy,  an  internal  state  vari¬ 
able  accounting  for  damage  or  defects,  and  a  logarithmic  measure 


j.D.  Clayton  /  International  Journal  of  Impact  Engineering  73  (2014)  56-65 


57 


of  thermoelastic  strain  (specifically,  the  logarithm  of  the  right 
thermoelastic  stretch  tensor). 

A  few  prior  works  have  used  the  multiplicative  decomposition 
to  describe  fracture  processes  in  brittle  solids  [14,15],  but  these 
focused  on  free  energy  rather  than  internal  energy,  and  used  Green 
elastic  strain  rather  than  logarithmic  strain.  The  latter  strain  proves 
more  convenient  for  analysis  of  planar  shocks  in  isotropic  solids, 
which  is  undertaken  in  the  second  part  of  this  paper  (§3).  Analysis 
here  follows  previous  works  dealing  with  elastic-plastic  metals 
[16,17],  the  primary  difference  being  that  in  brittle  solids,  pressure 
dependent  yield  strength  is  considered,  which  complicates  the 
solution.  In  this  work,  it  is  shown  how  the  Rankine-Hugoniot 
conditions  and  the  constitutive  equations  for  pressure  and  shear 
strength  can  be  reduced  to  a  set  of  coupled  algebraic  equations  that 
can  be  solved  for  the  material  state  downstream  of  the  shock,  given 
the  upstream  state  and  one  input  variable  describing  the  down¬ 
stream  state  (e.g.,  specific  volume  or  shock  pressure).  Pre-stress  in 
the  material,  prior  to  shock  compression,  is  also  newly  considered 
in  the  present  analysis. 

In  the  third  part  of  this  paper  (§4),  model  parameters  are 
populated  for  titanium  diboride  (TiB2)  and  solutions  are  generated 
for  the  shock  problem,  to  axial  shock  stresses  on  the  order  of 
50  GPa.  Titanium  diboride  is  a  strong  ceramic  whose  grains  have 
hexagonal  crystal  structure.  Under  impact  loading,  fracture  is 
thought  to  be  the  dominant  source  of  inelasticity  [4,6,7],  though 
pore  collapse  [5]  and  dislocation  motion  [18,19]  may  take  place  at 
high  pressures.  Multiple  investigations  have  discovered  that  tita¬ 
nium  diboride  demonstrates  a  rather  unique  “double  yield”,  which 
can  be  inferred  from  inflection  points  in  the  particle  velocity  history 
[5,20-22].  Values  of  the  HEL,  and  associated  physical  mechanisms, 
can  vary  substantially  with  microstructure  (grain  size,  porosity, 
impurities,  etc.  [5]),  which  in  turn  depends  on  processing  route. 
Spall  strength  degrades  rapidly  in  the  material  when  shocked 
above  the  HEL,  providing  indirect  evidence  that  fracture  is  the 
source  of  initial  yield  [4,7]. 

Key  new  contributions  of  the  present  paper,  in  the  context  of 
prior  work,  are  summarized  as  follows.  The  constitutive  model  and 
analytical  solution  consider  brittle  ceramics  with  pressure  depen¬ 
dent  strength,  in  contrast  to  [16,17]  focused  on  ductile  metals  with 
pressure  independent  yield.  In  contrast  to  [23],  which  modeled 
anisotropic  single  crystals  using  logarithmic  theory,  the  present 
work  considers  isotropic  polycrystals  (for  which  some  algebra 
simplifies  considerably)  with  possible  lateral  pre-stress,  with  the 
latter  condition  not  analyzed  previously.  Eurthermore,  the  current 
particular  application  to  titanium  diboride,  which  demonstrates 
the  aforementioned  unusual  double  yield  mechanism,  is  also  a  new 
contribution. 

Notation  of  continuum  physics  is  used:  vectors  and  tensors  are 
generally  written  in  bold  italic,  with  (•)^,  (•)“\  and  denoting 
transposition,  inversion,  and  the  inverse-transpose.  Individual 
components  of  vectors  and  tensors  are  referred  to  a  Cartesian  frame 
and  are  written  in  plain  italic,  with  summation  over  repeated 
indices.  The  scalar  product  of  two  second-order  tensors  is 
A:B  =  AijBij. 

2.  Finite  deformation  theory 

General  theory  is  developed  in  §2.1.  This  theory  is  specialized  to 
isotropic  thermoelastic  response  in  §2.2. 

2.1.  General  theory 

Let  X  =  x{X,t)  denote  spatial  coordinates,  at  time  t,  of  a  material 
particle  with  reference  coordinates  X.  Assuming  x  is  differentiable, 
the  deformation  gradient  is 


F  =  VoX  =  F^F°,  (2.1) 

where  Vq  denotes  the  referential  gradient  (e.g.,  in  Cartesian  co¬ 
ordinates,  Fij  =  djXi  =  dXijdXj)  and  and  denote  deformation 
“gradient”  mappings  associated  with  thermoelasticity  and  defects 
(e.g.,  cracks,  dislocations,  etc.),  though  neither  of  the  latter  gener¬ 
ally  anholonomic  mappings  [24]  need  be  compatible  (i.e.,  a  true 
gradient  of  a  vector  held).  Each  is,  however,  presumed  to  have 
positive  determinant.  In  the  absence  of  discontinuities,  body  forces, 
and  heat  sources/conduction,  the  usual  local  balance  laws  of  con¬ 
tinuum  mechanics  are  [8] 

Po  =  pJ,  V(7  =  pv,  (7  =  (7^,  U  =  (7:Vv,  podv>0.  (2.2) 

Reference  and  spatial  mass  densities  are  po  and  p,  volume  ratio  is 
J  =  detF,  Cauchy  stress  is  (7,  particle  velocity  is  u  =  x,  internal  energy 
per  unit  reference  volume  is  U,  absolute  temperature  is  6,  and  en¬ 
tropy  per  unit  reference  volume  is  r\.  The  spatial  gradient  is  V  [i.e., 
V/;(  )  =  0fc(  )  =  a(  )/0X/;],  and  time  derivatives  are  taken  at  flxedX. 
Helmholtz  free  energy  density  isW^U-dij. 

Applying  the  polar  decomposition  to  thermoelastic 
deformation, 

F^  =  =  F^^F^  =  =  R^^.  (2.3) 

Several  thermoelastic  strain  measures  are  deflned  for  later  use: 

=  =  Inl/^  =llnC^,  e^  =  lnV^.  (2.4) 

Definitions  and  identities  for  the  logarithm  of  a  tensor  are  dis¬ 
cussed  in  Ref.  [25].  Thermoelastic  volume  change  is 

f  =  detF^  =  detU^  =  detV^  =  det|^2F^  ^  ^  ^  exp(^tre^^ 

=  exp  • 

(2.5) 

Denote  by  f  a  generic  internal  state  variable  associated  with 
evolution  of  microstructure,  e.g.,  accumulated  cracks,  voids,  or 
dislocations  in  the  material.  Here,  f  is  assumed  a  scalar,  but 
generalization  to  higher-order  tensors  and/or  multiple  state  vari¬ 
ables  poses  no  difficulties.  Assuming  uniform  properties  in  the 
reference  state,  internal  energy  density  is  of  the  general  forms 

U  =  u(f\  v,  f)  =  U(£^  V,  f)  =  u{e^,  V,  f)  •  (2.6) 

Eorm  U,  incorporating  elastic  Green  strain  F^,  is  the  traditional 
choice  for  elastic  and  elastic-plastic  crystals  [8,9,26,27].  Eorm  U 
has  been  rarely  used  for  anisotropic  solids,  two  exceptions  being 
the  analysis  of  higher-order  elastic  moduli  in  Ref.  [28]  and  a  recent 
study  of  shock  compression  of  sapphire,  diamond,  and  quartz  single 
crystals  in  Ref.  [23].  Deflne  thermodynamic  conjugate  forces  as 

S  =  dU/dE^,  S  =  dU/de^,  ?  =  -aU/8f  =  -8i//8?.  (2.7) 

Using  (2.2)  and  considering  admissible  thermomechanical 
processes,  the  following  constitutive  laws  can  be  derived  consis¬ 
tently  with  the  first  and  second  laws  of  thermodynamics  (see  e.g. 
Ref.  [8]): 

^  =  jE-1pEsF^T  ^  jE-\pE  ET^  g  ^  qJJ ^  qQ 

ja  :  (F^f^F-1)  +  Cl  >  0. 

(2.8) 
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In  the  hyperelastic  law  in  the  first  of  (2.8),  the  following  fourth- 
order  tensor  is  used  [25]: 


M  = 


8lnC^ 

8C^ 


1  —  1  i  — 1  1— 1.W 


In  Af  —  In  Aj 
Af  -  Af 


PiSPf. 


(2.9) 


Here  Af  =  (xf)^  are  the  principal  values  of  C^, 
and 


Pi=  n  (2-10) 

Noting  that  principal  stretches  Af  are  eigenvalues  of  (and  V^), 
c^  =  ^P,lnAf  =  l^P,ln/lf.  (2.11) 

i=l  i=l 

Defining  specific  heat  at  constant  strain  as  c  =  dU/df)  =  -Sd^Wj 
d8^,  the  balance  of  energy  becomes 


8ora 


(9  877  Y  _  ( 


and 


Uo  =  U{0,'no),  Ca 


(2.19) 


(2.20) 


Oa(i 


JAl\  I  r  -  (  I 

<^4 J  ~  V0e£8e£8e£  J  'o 


(2.21) 


/O  dv\,  _  /d0\,  _  f  V 


(2.22) 


Letting  Cq  =  (dU/dO)!^  denote  a  constant  specific  heat  for  the 
unstrained  material,  the  strictly  entropic  contribution  to  internal 
energy  is  [29] 


Cd=J<T  :  (^F^f^F-i)  +  e(ds/de')  :  +  [?  -  ^(8?/8^)]| 

=  Ja  :  (F^iT’p-''')  +0^ds/d0^  :  ^  +  [C  -  6»(8C/66»)]f. 

(2.12) 

Generic  kinetic  equations  for  inelasticity  are  of  the  state 
dependent  form 

=  F^(fCv,?),  ?  =  ?(F^r,,Y.  (2.13) 

Let  Uq  denote  internal  energy  in  the  elastically  unstrained 
reference  state  defined  by  {E^,7])  =  (0,770)*  =  (0,770)*  and  en¬ 

tropy  change  from  this  reference  state  is  A77  =  77  -  770.  Temperature 
change  from  this  reference  state  is  -  Oq.  Let  Greek  subscripts 
denote  Voigt  notation  for  symmetric  indices,  e.g.,  {-)jj  ^ 


h  =  Co[exp(A77/co)  -  1]  =  +  +  +  •••■ 

(2.23) 

It  can  be  shown  [17,28,29]  that  second-order  isentropic  elastic 
constants  and  Griineisen  constants  Fa  should  be  equal  when  the 
reference  state  is  unstressed  for  £-based  theory  and  e-based  theory: 

Ca0  =  Cap  =  Cap,  Ta  =  Ta=  Fa-  (2.24) 

This  result  is  consistent  with  the  requirement  that  U^U  when 
strains  are  small.  Third-order  isentropic  constants  are  related,  in 
full  tensor  notation,  by  Ref.  [28] 

ClJKLMN  =  CjjklMN  +  2  (JijkLPqCpqmN  FJklMNPqCpqjj  FJmNIJPqCpqkL^  , 

(2.25) 


11^1,  22^2,  33^3,  23  =  32^4,  13  =  31^5, 

12  =  21^6. 

(2.14) 

The  following  Taylor  series  expansions  of  internal  energy  are 
used: 

U=Uo  +  CX  +  -  So  -Hv)] 

+/(f), 

(2.15) 


JlJKLMN  =  g  {^IK^jMhN  +  +  ^IL^JM^KN  +  ^IL^JN^KM 

+  ^IM^JKhN  +  ^IM^JL^KN  +  ^IN^JK^M  +  • 

(2.26) 

Henceforward  this  work  will  use  e-based  theory  and  Lf,  with 
(2.25)  needed  only  to  convert  experimentally  reported  values  of 

Capj  to  CaPj‘ 

2.2.  Isotropic  thermoelasticity 


aE  i  ^  r  oEaE  ^ 


U=Uo  +  C„e^  -  8o 


Tael^v-h{v) 


+/(?)• 


(2.16) 


Thermoelastic  relations  simplify  considerably  when  the  solid  is 
elastically  isotropic.  Hyperelastic  law  (2.8)  reduces  to  the  following 
relation  between  principal  components  of  Cauchy  stress  (Ti  and 
elastic  principal  stretch  components  Af  [17,30]: 


Material  coefficients  evaluated  at  the  unstressed  reference  state 
are 


Uo  =  U{0,vo),  Q  =  (8D/6£Y|o  =  0; 


(2.17) 


C  -  (  )l  C  -(  )l 


(2.18) 


_  1  EdU  _  1  6L/  _  1  8D 

'"'■“4  ’8Af 


1,2,3). 


(2.27) 


Note  from  (2.11)  that  ef  are  principal  values  of  e^,  equal  to  log¬ 
arithmic  principal  stretches  In  Af .  Isotropic  second-  and  third-order 
elastic  constants  are  of  the  following  forms  (C  or  C): 


ClJKL  =  ^^IJ^KL  +  FipIK^JL  +  ^IE^Jk)  , 


(2.28) 
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QjKLMN  =  I^IJ^KL^Mn]  +  ^2  l^Iji^KM^LN  +  ^KN^Lm)  +  ^Kl{^IM^JN 

+  %^;m)  +  ^mn{^ik^jl  +  ^il^jk)]  +  ^3  [^ik{^jm^ln 
+  ^JNhlVl)  +  ^jd^IM^KN  +  %%m)  +  ^Ili^JM^KN 
+  ^JN^Km)  +  ^JKi^IMhN  +  %^Lm)]  • 

(2.29) 

Second-order  constants  obey 
Cii=A  +  2^,  C|2  =  C44  =  V  =  Xl{2X 2iii), 

K  =  A  + 

(2.30) 

Shear  modulus  is  bulk  modulus  is  /C,  and  Poisson's  ratio  is  j'. 
Third-order  constants  obey  the  following  relations  or 

Qii  =  +  6za2  +  Qi2  =  +  2i^2,  Ci23=^i, 

Ci44  =  iA2,  C155  =  i^2  +  2za3,  C456  =  Z^3. 

For  isotropic  materials,  (2.25)  reduces  to 

3 

=  h=^3^2^'  (2.32) 

For  hydrostatic  elastic  loading  {a  =  -pi,  F  =  pressure 

derivatives  of  tangent  bulk  modulus  B  and  tangent  shear  modulus  G 
at  the  reference  state  are  [29,31,32] 

B'o  =  {dB/dp)\^^^  =  -1  (f,  +  2^2  +  |i73^ 

=  -J(^i,,+2u2+lh^+2.  (2.33) 

g;,  =  (dG/dp)|^^^  =  -1  (^2  +  +  5-“)  - 1 

= (^^2 +|»'3  -  A -  1.  (2.34) 

The  above  relations  hold  for  isothermal  or  isentropic  elastic 
constants,  but  the  notation  refers  herein  to  isentropic  constants, 
consistent  with  definitions  in  §2.1.  The  Gruneisen  tensor  is  spherical 
for  isotropic  solids  [8]: 

Fjj  =  rdjj  =  (3aK/Cp)djj.  (2.35) 


moderately  large  strains  has  been  demonstrated  elsewhere  [33]. 
Finally,  logarithmic  equation  of  state  (2.36)-which  is  derived 
naturally  from  the  complete  logarithmic  nonlinear  theory-has 
been  shown  [34]  to  be  highly  realistic  in  geophysics  applications 
spanning  small  to  very  large  pressures.  Justification  of  use  of  the 
logarithmic  theory  (e.g.,  in  contrast  to  linear  elasticity)  to  describe 
titanium  diboride  in  particular  will  become  clear  in  §4  in  the 
context  of  comparisons  of  constitutive  model  predictions  to 
experimental  compression  data. 

3.  Shock  compression  analysis 

Governing  equations  for  planar  shock  loading,  i.e.,  jump  condi¬ 
tions,  are  reviewed  in  §3.1.  Application  of  the  present  finite  strain 
constitutive  model  is  presented  in  §3.2,  and  the  method  of 
analytical  solution  is  outlined  in  §3.3. 

3.1.  Governing  equations 

Consider  a  continuous  cylinder  of  material  through  which  a 
planar  shock  moves,  in  the  Xi -direction,  with  natural  velocity  Let 
superscripts  +  and  -  label  quantities  in  the  material  ahead  (i.e., 
upstream)  and  behind  (i.e.,  downstream)  from  the  shock.  Let  [(•)]] 
and  ((•))  denote  the  jump  and  average  of  a  quantity  across  the 
shock: 

I(-)l  =  (•)■-(•)+.  <(•))  =  ^[(•)"  +  (-)^].  (3.1) 

Let  n  be  a  unit  normal  vector  to  the  planar  shock,  i.e.,  n  =  dx/dxp 
The  only  nonvanishing  component  of  particle  velocity  is  v  =  v-n. 
The  Cauchy  stress  component  normal  to  the  shock  front  is 
(T  =  <7  :  (ii(8)ii)  =  (Tn.  The  relative  velocity  of  the  material  with 
respect  to  the  shock  is  v  =  u  -  Let  u  =  U/po  denote  internal 
energy  per  unit  mass.  Appropriate  forms  of  the  Rankine-FIugoniot 
conditions  for  conservation  of  mass,  momentum,  and  energy  are 
[16] 

M  =  0,  (3.2) 

M-p4i^]  =  0,  (3.3) 

lpv(u  +  -  av]i  =  0.  (3.4) 


Here,  a  is  the  linear  coefficient  of  thermal  expansion,  and  spe¬ 
cific  heat  at  constant  pressure  in  the  reference  state  is 
cp  =  co(l  +  BaTdo).  Under  hydrostatic  elastic  and  isothermal 
loading,  applying  (2.27)  and  (2.33),  the  present  model  degenerates 
to  the  following  pressure-volume  equation  of  state: 


P 


-dW  /dj 


l-l(fij,-2)ln; 


(2.36) 


where  Kq  =  /Cco/cp  is  the  isothermal  bulk  modulus. 

General  benefits  of  the  logarithmic  finite  strain  theory  can  be 
explained  as  follows.  The  principal  Cauchy  stress-logarithmic 
elastic  strain  relations  (2.27)  demonstrate  a  particularly  simple 
and  convenient  form  for  analysis  of  isotropic  materials  [33].  The 
logarithmic  elastic  and  plastic  strains  resulting  from  a  multiplica¬ 
tive  decomposition  of  the  deformation  gradient  (2.1)  combine  to  a 
basic  additive  split  of  the  total  logarithmic  strain  tensor  in  (3.10)  of 
§3.2,  which  facilitates  mathematical  analysis  of  the  uniaxial  shock 
problem  [16].  Accuracy  of  the  elastic  theory  for  describing  the 
response  of  a  wide  class  of  different  materials  subjected  to  various 
loading  protocols  (extension,  compression,  and  shear)  up  to 


The  material  need  not  be  deformed  uniaxially  according  to  these 
conditions,  but  the  shock  velocity  and  particle  velocity  must  both 
be  rectilinear  in  the  xi -direction  so  that  only  normal  traction  is 
discontinuous.  Therefore,  these  equations  can  apply  for  shock(s) 
passing  through  a  pre-stressed  material,  as  will  be  considered  later. 
Adiabatic  conditions  have  been  assumed  [16]:  heat  conduction  is 
not  included,  leading  to  entropy  production  requirement 

M  >  0.  (3.5) 

The  shock  process  is  neither  isentropic  nor  isothermal,  in  gen¬ 
eral.  Using  (3.2)  and  (3.3),  energy  conservation  condition  (3.4)  can 
be  rewritten  as  [16] 

M  =  {ct)U/pI  (3.6) 

Assume  that  the  upstream  state  and  shock  velocity  ^  are 
known.  The  downstream  state  can  be  defined  by  the  set  of  variables 
(u“,p“,(7“,u“).  The  Rankine-Hugoniot  conditions  provide  three 
equations  for  determining  this  state;  in  order  to  fully  determine  the 
downstream  state,  a  fourth  equation  must  be  supplied  by  the 
constitutive  model  for  the  shocked  material,  or  one  more  of  the 
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downstream  variables  must  be  known,  e.g.,  observed  from  exper¬ 
iment  (typically  the  particle  velocity  v~  in  addition  to  ^). 


3.2.  Shock  compression  analysis  of  an  isotropic  brittle  solid 

Consider  a  brittle  material  whose  thermoelastic  response  is 
described  by  the  theory  in  §2.2.  Assume  that  the  material  may  be 
statically,  isothermally,  and  elastically  pre-deformed  prior  to  shock 
compression,  with  the  following  initial  elastic  deformation 
gradient  in  3x3  matrix  form: 


U  =  u(ef,  v)  =^Cijefef  +  Ic^efe/e^  -  doFv  \nj 

+  0o(^V  +  j ^  (^2) 

+  2X6^62  +  (^0^1+  ^2  +  3  ^3^  (cf  j 


^  (^2)  +  (^1  + 


4_  .  8_ 

j!'!  +4i'2  +2!'3 


^1 

0 

0  ■ 

(3.7) 

+  ^0  (  ^  +  /^O  ) 

0 

^2 

0 

V  2  /  7 

0 

0 

^2_ 

(3.14) 


This  applies,  for  example,  to  a  cylindrical  specimen  under  uni¬ 
form  lateral  pressure.  Subsequently,  let  the  specimen  be  subjected 
to  uniaxial  strain  loading  of  the  form  Fi,  which  in  general  may 
consist  of  ramp  loading  and/or  multiple  planar  shocks,  and  may  be 
elastic  and  inelastic: 


Fi 


Ai  0  0 

0  1  0 
0  0  1 


The  total  deformation  gradient  of  (2.1 )  is  then 


(3.8) 


The  following  assumptions  have  been  made.  First,  internal  en¬ 
ergy  increase  associated  with  generation  of  new  crack  surfaces  is 
assumed  to  offset  internal  energy  decrease  associated  with  reduc¬ 
tion  in  elastic  moduli  as  damage  progresses,  the  latter  phenomenon 
being  omitted.  Therefore,  dependence  of  U  on  f  is  omitted  {/=  0). 
Second,  in  the  unstressed  reference  state,  conditions  Uq  ^  0  and 
7]o  ^  0  are  chosen,  an  assumption  which  does  not  affect  the  me¬ 
chanical  response  that  depends  only  on  changes  in  U  and  p  and  not 
their  absolute  values.  Third,  entropic  contribution  (2.23)  is  trun¬ 
cated  at  second  order.  From  (2.27),  principal  Cauchy  stress  com¬ 
ponents  are 


AiAi 

0  0  ■ 

FiFo  = 

0 

42  0 

=  F^ 

0 

0  A2 

0 

0  ' 

‘Af  0 

0  ' 

0  Af 

0 

0 

0 

0  0 

^2  J 

0  0 

^2_ 

(3.9) 


0-1  =J 


(A  -h  2^)0|  -h  27^2  +  [  2  “1“  ^^2 


(3.15) 


Since  F  is  symmetric,  F  ^  U  (and  similarly  for  its  elastic  and 
inelastic  parts),  so  the  total  logarithmic  strain  reduces  to  the  ad¬ 
ditive  form 

\nU  =  InF  =  InF^  +  InF^  =  +  e^.  (3.10) 

Flenceforth  it  is  assumed  that  inelastic  deformation  is  isochoric, 
corresponding  primarily  to  mode  11/111  opening  and  sliding  of 
micro-cracks.  Inelastic  volume  changes  associated  with  pore 
collapse  and  crack  opening/dilatation  are  assumed  to  be  offsetting. 
Thus,;  =  Af  (Af)2  and 

tre°  =  ln(detF'^)  =  0=>ef  =  -^e?,  (3.11) 

where  ef  =  In  Af .  Assume  that  Ai  and  A2  are  known  a  priori,  e.g., 
from  a  static  elasticity  calculation.  Then  (3.9)— (3.11)  can  be  com¬ 
bined  to  give  the  following  expressions  for  logarithmic  elastic  strain 
components  in  terms  of  e\  =  InAi  and  a  single  inelastic  strain 
component  e^: 

ef  (ei ,  ef )  =  ei  -h  In  /li  -  ef ,  ef  (ei .  ef )  =  In  /I2  +  ^ef  • 

(3.12) 

Total  volume  change  from  the  undeformed,  unstressed  state  is 
ln/  =  ln;  =  ei -Hln/\i -H21n/l2.  (3.13) 


(72  =  (73 

=  7  ^  I^Ae^  -h  (2A  -h  'lp)e\  -\-  +  ^2^  (^0  ^^2 

+  ‘4i>3^^e2^  (2i>i  2i>2^C^C2  —  OoFri  . 

(3.16) 

Cauchy  pressure  is 


P 


1  2 
“  3^2 


14.  8. 

+  (  +^^^2  +  3^^3 


^^2)  +  (5^1  +|^2+|^'3)(ef) 

^(^2)  +  -i-|i>2^efef  -  6»ol’’7 


(3.17) 


Let  (7-  =  (7j  +  p  denote  deviatoric  stress  components.  The 
following  quantities  are  associated  with  shear  stress: 


ff2) 


-|(<^i+p)  =  ~<^'i  =\o'2=\7‘hy^ 


(3.18) 


The  thermoelastic  response  is  now  addressed.  Internal  energy 
(2.16)  becomes,  for  the  present  situation  where  all  strain  compo¬ 
nents  are  principal  strains  (ij,/<  =  1,2,3): 


where 72  =  1  /2[((7'^)^  +  (cr^)^  +  ((73)^],  and  Oe  is  the  effective  (Mises) 
stress,  equal  to  the  applied  stress  in  a  static  uniaxial  stress  exper¬ 
iment.  From  (3.15)  and  (3.16), 
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r  =  -J 


-1 


Ailef-ef 


+  (v2  +  2V3)  -  (2^2  +  2vs')  (^efj 


F  F 


(3.19) 


closed  form,  but  they  can  be  evaluated  exactly  in  straighforward 
manner  via  simple  iteration  using  a  computer  program.  Once  shear 
stress  and  pressure  are  computed,  shock  velocity  ©  and  down¬ 
stream  particle  velocity  v~  can  be  obtained  from  the  Hugoniot 
equations  for  mass  and  momentum  conservation,  (3.2)  and  (3.3), 
leading  to 


In  the  present  application  of  the  general  theory  of  §2.1  to  brittle 
isotropic  solids  under  planar  shock  compression,  a  rate  indepen¬ 
dent  constitutive  relation  is  used  to  specify  inelastic  deformation  in 
lieu  of  (2.13).  A  general  yield  criterion  is 

r  =  lY;  y  =  y(c£.„.f),  f  =  f(e?).  (3.20) 

Dynamic  yield  strength  Y  can  depend  strongly  on  pressure 
as  mentioned  earlier  an  important  feature  of  brittle  ce¬ 
ramics  that  does  not  typically  apply  to  ductile  metals  [16,17].  In¬ 
ternal  variable  f  represents  damage  in  the  material  and  is  assumed 
to  increase  monotonically  with  plastic  strain  until  saturation  (i.e.,  a 
percolation  limit  of  micro-cracks,  or  pulverization).  More  general 
relations  would  be  needed  to  address  non-monotonic  loading 
conditions.  As  a  more  specific  example,  the  following  strength 
expression  Y  =  Y(p^)  is  considered: 

y  =  yo[i  +  x(p,/yo-PH/yo)'"],  p,(e^)  =P-6orvfj- 

(3.21) 

Here,  explicit  dependence  on  f  is  excluded,  and  strength  de¬ 
pends  on  isentropic  pressure  which  eliminates  coupling  with 
entropy  to  enable  an  exact  analytical  solution  to  the  shock  problem, 
as  discussed  later.  Physically,  this  simplification  corresponds  to 
offsetting  entropic  pressure-shear  hardening  and  thermal  soft¬ 
ening  associated  with  defect  kinetics.  Yield  stress  at  the  HEL  is 
Yo  >  0,  isentropic  pressure  at  the  HEL  is  ph,  and  x  ^nd  m  are  pa¬ 
rameters.  Eor  example,  other  models  of  pressure  dependent 
strength  of  brittle  materials  essentially  prescribed  m  =  1/2  [2]  and 
m  =  1  [12] . 

3.3.  Analytical  solution 

Independent  solution  variables  are  the  set  (p,T,  6^,77,61),  with 
ei  =  InAi  prescribed  as  the  loading  parameter.  If  the  material  is 
subjected  to  static  lateral  pre-stress  of  magnitude  ao  =  -a2  =  -^^3. 
constants  A\  and  A2  are  determined  from  simultaneous  solution  of 
(3.15)  and  (3.16),  with  =  0,  ef  ^InAi,  ef  ^lnA2,  p  ^  0,  and  the 
isothermal  modulus  X  ^  Xq  =  Kq  -  213 p.  Elastic  strains  can  then  be 
written  as  functions  of  e\  and  using  (3.12).  Equations  (3.17)  and 
(3.19)  can  be  written  in  the  form  p  =  p(ei ,  p)  and  r  =  T(ei ,  e^), 
providing  two  independent  equations.  Combining  (3.20)  and  (3.21 ) 
gives  a  third  independent  equation  of  the  form  r  =  1  /2Y(ei ,  e^).  A 
fourth  independent  equation  is  the  energy  conservation  law  across 
the  shock,  (3.6),  which  can  be  written 

m  =  {<Ttm  =  -{^r+p)M.  (3.22) 

where  internal  energy  in  (3.14)  can  be  expressed  in  the  form 
f;  =  f;(ei,e?,p). 

The  analytical  solution  proceeds  as  follows.  Equating  r{e-i ,  e^)  = 
1  /2Y(ei ,  e^)  gives  a  quadratic  algebraic  equation  that  can  be  solved 
for  =  e^(ei ),  and  then  r[e-i ,  e^(ei )]  and  elastic  strain  components 
can  all  be  calculated  explicitly  in  terms  of  loading  variable  e\. 
Equations  (3.17)  and  (3.22)  can  then  be  solved  simultaneously  for 
p(ei)  and  p(ei),  assuming  that  the  upstream  state  is  known  a  priori 
in  (3.22).  The  resulting  solutions  are  too  lengthy  to  express  in 


®  =  {  ((7^  -  (7+ j  /  [po  (■^r  -  ■^1")]  }  ^  ®  (a^  - 

(3.23) 

Here,  po  is  the  mass  density  after  possible  pre-stress  but  prior  to 
shock  compression.  Eor  a  single  plastic  wave,  typically  it  is  assumed 
that  the  upstream  state  is  elastic  [i.e.,  =  0]  but  has  been 

stressed  to  the  HEL  by  the  elastic  precursor.  The  upstream  state  is 
then  the  set  of  variables  (p+,T+,r7+,e|),  where  the  axial  positive 
compressive  stress  at  the  HEL  is  (Jh  =  p”^  +  4/37“^  and  the  corre¬ 
sponding  strain  induced  by  the  shock  is  =  In  •  The  shock 

velocity  and  particle  velocity  of  the  elastic  precursor  are 

®£  =  -  exp(eH))]}’^^  i-e  =  ®£[1  -  exp(eH)]. 

(3.24) 

If  Specific  volume  at  the  HEL  (i.e.,  en)  is  known  for  an  initially 
unstressed  material,  Yq  can  be  determined  from  T'^(eH)  =  1/ 
2Y+(eH)  =  Yq.  Effects  of  entropy  production  are  typically  small  even 
for  moderate  elastic  shocks  [27,29]  and  can  thus  usually  be  omitted 
for  the  precursor,  i.e.,  ^  po  =  0-  Temperature  is 

0  =  dU /dn  =  0o[\-  r{e\  +  2e|)  +  v/Cq] •  (3.25) 

The  constitutive  model  framework  of  §2  and  method  of  solution 
to  the  shock  problem  outlined  above  in  §3,  while  applied  specif¬ 
ically  to  titanium  diboride  in  §4,  are  general  enough  to  be  poten¬ 
tially  applicable  to  a  large  number  of  solids  undergoing  finite 
deformation,  so  long  as  their  response  can  be  treated  as  isotropic 
hyperelastic-plastic,  with  shear  strength  possibly  depending 
strongly  on  pressure  in  the  inelastic  regime.  Eor  example,  the  shock 
response  of  brittle  polycrystalline  rocks  and  minerals  [35],  as  well 
as  other  ceramics  such  as  silicon  carbide  [36]  that  yield  predomi¬ 
nantly  by  micro-cracking,  may  be  considered.  The  analysis  could 
also  be  applied  towards  description  of  the  shock  response  of  metals 
whose  inelasticity  is  governed  by  dislocation  slip,  as  considered 
previously  in  Ref.  [17]  and  for  which  pressure  dependence  of 
strength  may  be  reduced  or  omitted.  However,  more  sophisticated 
numerical  techniques  such  as  the  plane  wave  method  [37]  are 
needed  to  address  effects  of  dislocation  interactions  and  slip  rates 
on  shock  profiles. 


4.  Shock  response  of  titanium  diboride 

Properties  of  titanium  diboride  are  described  in  §4.1.  Solutions  to 
the  shock  problem  are  presented  and  discussed  in  §4.2. 


4.1.  Properties 

Titanium  diboride  (TiB2)  is  a  hard  crystalline  ceramic  whose 
grains  have  hexagonal  symmetry.  Average  grain  sizes  in  polycrystals 
have  been  reported  to  vary  from  approximately  5  to  50  ^m,  and 
initial  densities  from  -95%  to  nearly  100%  of  theoretical  density  of 
4.52  g/cm^  [5,38].  Properties  needed  for  application  of  the  analysis 
and  solution  of  §3  are  listed  in  Table  1.  Ambient  density,  second-order 
elastic  constants,  and  thermoelastic  constants  are  obtained  directly 
from  the  literature.  Poisson's  ratio  of  this  particular  material 
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(corresponding  to  that  manufactured  by  Cercom  [4,5])  is  =  0.049, 
somewhat  lower  than  that  of  other  samples  reported  in  Ref.  [4]. 

Direct  experimental  measurements  of  third-order  elastic  con¬ 
stants  of  TiB2  have  apparently  not  been  reported.  However,  values 
of  Bq,  Gq,  and  C[  have  been  measured  ultrasonically  [38],  where  the 
first  two  are  described  by  (2.33)  and  (2.34)  and  C[  is  the  pressure 
derivative  of  the  tangent  longitudinal  elastic  modulus  in  the 
reference  state.  The  latter  can  be  related  to  second-  and  third-order 
elastic  constants  as  [31,39] 


C[  =  -  (si'i  +  lOj/2  +  81-3  -  lOA  -  Ufj.J/(3K)  -  (A  +  2fj.) 

/(3K)--l. 

(4.1) 

Given  Gq,  C[,  and  the  second-order  elastic  constants.  Equa¬ 
tions  (2.33),  (2.34),  and  (4.1 )  cannot  be  solved  simultaneously  for  all 
three  independent  third-order  constants  since  these  equations  are 
not  linearly  independent.  Instead,  herein  the  experimentally  re¬ 
ported  axial  stress  at  the  HEL,  (Jh,  for  an  initially  unstressed 
isotropic  sample. 


-  Jh’ 


Cll^H  +  2^iii^h 


-1 


=  -Jh 


(2  +  2^)eH  ■ 


-Ai's  ]e^ 


2^  +3p2 


(4.2) 


with  experimentally  reported  density  ratio  po/p  =  Jh  =  exp(eH)  is 
used  with  (2.33)  and  (2.34)  to  yield  the  third-order  constants 
shown  in  Table  1.  Explicitly,  equations  for  third-order  constants  are 

v,=6K{G'o  +  -i)-2{3X  +  5^i)  +  C,u,  (4.3) 


h=j2  -  31<:  (3BJ)  +  1 6G^,  +  1 0)  +  48 A  +  80^  -  9Ci  1 1 
P3=T  3/C(B[,  +  4G[,  +  2)-4(3A  +  5m)  +  3Ci„  , 


(4.4) 

(4.5) 


Fig.  1.  Axial  (Hugoniot)  stress  P  =  -a-i\  model  and  experiment  [5].  Static  lateral  pre¬ 
stress  is  a  o',  volume  after  pre-stress  but  before  shock  compression  is  Vq. 


idealization  has  been  shown  acceptable  in  previous  analysis  of 
elastic  shocks  up  to  the  HEL  in  quartz,  sapphire,  and  diamond  [29]. 

As  noted  already,  titanium  diboride  exhibits  an  inelastic  wave 
structure  following  the  elastic  precursor  that  can  be  associated 
with  two  “HEL”  states  [5,20—22],  herein  labeled  HI  and  H2.  The 
second  is  not  truly  an  elastic  limit  since  inelasticity,  in  the  form  of 
micro-cracking,  commences  in  the  material  after  initial  yield  [7,21  ]. 
A  failure  wave  associated  with  lateral  stress  degradation  has  also 
been  reported  at  shock  pressures  between  first  and  second  yield 
[21].  The  presence  of  two  yield  points  can  be  inferred  from  in¬ 
flections  in  the  particle  velocity  history  and  also  is  apparent  from 
some  experimental  Hugoniot  pressure-volume  curves  [5,20].  To 
capture  the  effect  in  the  model,  the  pressure  dependent  strength 
equation  is  applied  separately  to  regimes  shocked  to  Hugoniot 
states  below  and  above  the  second  yield  point,  with  different 
constants  used  in  (3.21 )  for  each  regime.  Let  P  =  -an  =  denote 
the  axial  stress  (i.e.,  shock  pressure),  positive  in  compression,  along 
the  Hugoniot.  Let  am  and  am  denote  values  of  P  at  the  first  and 
second  yield  points.  Eor  titanium  diboride,  yield  strength  of  (3.21 )  is 
prescribed  as 


Cm  =  -  2  [(A  +  2ix)/eu  +  an  exp(eH) / e^] ;  (4.6) 

where  au  and  en  are  taken  from  experiment  “SNL-2”  of  Grady  [5]  on 
Cercom  material;  ultrasonic  and  shock  data  on  this  material  are  also 
tabulated  in  Ref.  [4].  Equations  (4.2)  and  (4.6)  are  approximate 
because  entropy  rise  for  the  elastic  process  is  omitted.  This 


Table  1 

Physical  properties  of  TiB2  {do  =  295  K;  po  in  g/cm^;  Co  in  MPa/K;  K,  p,  (Th,  Yq  in  GPa). 


Label 

Description 

Value 

Reference 

Po 

ambient  mass  density 

4.51 

[4] 

K 

isentropic  bulk  modulus 

193 

[4] 

P 

elastic  shear  modulus 

249 

[4] 

r 

Griineisen  constant 

1.1 

[38] 

Co 

specific  heat 

2.76 

[40] 

B'o 

pressure  derivative 
of  bulk  modulus 

2.18 

[4] 

^0 

pressure  derivative 
of  shear  modulus 

2.53 

[38] 

O-H 

HEL  stress 

6.0 

[5] 

Ch 

HEL  strain 

-0.00965 

[5] 

(Yoh 

strength  2t  at  HEL 

6.16 

[5]  with  Equation  (4.7) 

(Yo)2 

strength  2t  at  second  yield 

10.96 

[5]  with  Equation  (4.7) 

Xi 

strength-pressure  constant 
below  second  yield 

1.0 

model  fit 

X2 

strength-pressure  constant 
above  second  yield 

0.5 

model  fit 

y=(yo)i{i+xi 

y=(yo)2{i+x2 


/  /  1V2'| 

Pp/ (^o)i -Phi/(^o)i  I  for 

Pp/(^o)2-PH2/(^o)2  }  for  P>0'h2- 


(4.7) 


Thus,  m  =  1/2  for  the  response  up  to  second  yield  (e.g.,  as  in  Ref. 
[2]),  and  m  =  1  for  the  response  beyond  second  yield  (e.g.,  as  in  Ref. 
[12]).  Constants  (Yq)!  and  (Yo)2  correspond  to  strength  at  HI  and 
H2;  the  former  is  constrained  by  (Yo)i  =  2T(eHi),  where  r  is  known 
from  the  nonlinear  elastic  solution  at  HI,  and  the  value  of  the  latter 
is  fixed  to  ensure  continuity  of  the  Hugoniot  stress  between  the  two 
regimes  at  H2.  Pressures  pni  and  ph2  are  calculated  during  the 
analysis  and  are  not  free  parameters.  Thus,  essentially  only  two 
parameters  are  considered  adjustable:  xi  and  X2-  These  are  pre¬ 
scribed  to  fit  the  Hugoniot  stress-volume  data  of  [5]  for  initially 
unstressed  TiB2,  as  shown  later  in  Eig.  1.  Other  results  presented 


Table  2 

Error  in  model  predictions  versus  experimental  data  [4,5]  (no  pre-stress). 


VIVo 

Error  in  P  [%] 

Error  in  p  [%] 

Error  in  r  [%] 

0.983 

+10 

-1 

+12 

0.971 

-6 

-0 

-29 

0.943 

+1 

-1 

+5 

0.922 

+4 

-1 

+22 

0.920 

+7 

-1 

+29 
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later  for  mean  stress,  shear  stress,  inelastic  deformation,  temper¬ 
ature,  and  entropy  in  §4.2  (Figs.  2-4)  are  all  predictions  rather  than 
model  fits. 

For  samples  with  compressive  pre-stress  (i.e.,  ao  >  0),  inelas¬ 
ticity  and  strength  are  treated  as  follows.  Define  the  internal  state 
variable  associated  with  cumulative  inelastic  mechanisms  (i.e., 
fracture  and  possible  slip)  as 

?=  ||e»||  =-^3726?,  (4.8) 

where  (3.11 )  has  been  used.  Initial  yield  in  the  pre-stressed  material 
is  defined  by  the  condition  r  =  Vfff,  where  Y  is  given  by  the  first  of 
(4.7)  and  pni  is  the  isentropic  pressure  at  Fll  for  the  material 
without  pre-stress.  Since  increases  with  compressive  pre-stress 
at  the  same  value  of  ei,  the  shock  stress  and  compressive  strain 
associated  with  the  shock  at  HI  will  increase  with  increasing  pre¬ 
stress,  as  shown  later  in  §4.2.  Let  fc  be  the  value  of  ?  predicted  by 
the  model  at  H2  in  the  material  without  pre-stress,  which  corre¬ 
sponds  to  accumulation  of  a  critical  amount  of  damage.  Second 
yield  in  the  material  with  pre-stress  is  assumed  to  occur  when 
?  >  ?c;  with  increasing  compression  beyond  this  point,  the  second 
of  (4.7)  is  used  for  strength  in  the  pre-stressed  material.  Thus,  no 
additional  fitting  parameters  are  used  to  model  the  material  with 
pre-stress,  meaning  results  shown  later  in  Fig.  5  are  fully  predictive. 

The  present  analysis  suggests  that  use  of  nonlinear  elasticity  is 
essential  for  accurate  stress  predictions  of  TiB2  under  shock 
compression.  For  the  initially  unstressed  material  at  e\  =  -0.00965, 
the  present  nonlinear  elastic  model  gives  P  =  6.00  GPa, 
r  =  3.08  GPa,  and  p  =  1.90  GPa.  In  contrast,  at  the  same  strain  and 


with  \i  and  p  from  Table  1,  linear  elasticity  predicts  P  =  5.06  GPa, 
r  =  2.40  GPa,  and  p  =  1.86  GPa.  Restricting  attention  now  to  the 
pressure-volume  response  (i.e.,  hydrostatic  compression  curve), 
nonlinear  theory  in  (2.36)  predicts  p  =  22.5  GPa  at  V/Vo  =  0.9,  a 
result  verified  later  in  §4.2  as  highly  accurate  upon  comparison  with 
experiment.  In  contrast,  linear  elasticity  under-predicts  the  pres¬ 
sure  as  p  =  /C(l  -])  =  19.3  GPa  at ;  =  V/Vq  =  0.9. 

4.2.  Results 

Model  results  for  axial  stress  P  are  compared  with  shock  data  of 
Grady  [5]  in  Fig.  1.  Experimental  conditions  correspond  to  no  pre- 


Fig.  2.  Model  and  experiment  [4]:  (a)  mean  stress  (pressure)  p  versus  volume,  (b)  shear  stress  r  versus  volume,  and  (c)  shear  stress  versus  mean  stress. 
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stress  (do  =  0).  Agreement  between  model  and  experiment  is 
considered  good,  with  differences  between  model  and  experiment 
on  the  order  of  10%  or  less  in  Table  2.  [Errors  in  Table  2  are 
computed  as  (predicted  value  minus  experimental  value )/( average 
of  predicted  and  experimental  values).]  Changes  in  slope  of  P  versus 
volume  VIVq  =  po/p  are  evident  among  the  elastic  regime  (below 
HI),  the  first  plastic  regime  (between  HI  and  H2),  and  the  second 
inelastic  regime  (above  H2).  Compressive  pre-stress  leads  to  an 
increase  in  shock  pressure  P  at  the  same  volume  ratio. 

Predicted  mean  stress  (i.e.,  Cauchy  pressure  p,  positive  in 
compression)  is  compared  with  experimental  data  [4]  in  Fig.  2(a). 
The  upper  three  curves  are  results  from  shock  compression. 
Equation  (3.17).  The  model  hydrostat  is  computed  via  (2.36),  and  is 
lower  than  the  shock  compression  curve  at  null  pre-stress  pri¬ 
marily  because  the  hydrostat  does  not  include  pressure  rise  due  to 
entropy  generation  present  in  the  shock  calculation.  Agreement 
between  model  and  experiment  is  excellent,  with  with  differences 
(i.e.,  relative  error)  on  the  order  of  1%  or  less  in  Table  2. 
Compressive  pre-stress  results  in  a  predicted  increase  in  mean 
stress  under  shock  compression.  It  has  been  suggested  [5]  that 
pore  collapse  may  be  an  important  source  of  inelasticity  in  tita¬ 
nium  diboride,  particularly  in  samples  of  lower  initial  density. 
Here,  material  of  very  high  nominal  density  is  considered  (99.8%  of 
theoretical  density),  justifying  omission  of  the  effect  of  porosity  on 
compressibility. 

Predicted  shear  stress  is  compared  with  experimental  data  [4]  in 
Fig.  2(b),  where  the  latter  was  determined  from  the  offset  between 
longitudinal  stress  and  the  hydrostat  [4].  Model  results  are  found 
via  (3.19).  Predicted  trends  are  in  reasonable  agreement  with  the 
data  of  [4],  with  predicted  values  somewhat  larger  than  the  test 
data  in  the  second  inelastic  regime  (above  H2),  as  is  clear  from 
Table  2.  The  model  predicts  small  overall  effects  of  pre-stress  on 
shear  stress,  but  r  does  increase  with  gq  at  large  compression  due  to 
increasing  p^. 

As  observed  in  Fig.  2(c),  shear  stress  r  increases  significantly, 
and  in  a  nonlinear  manner,  with  increasing  compressive  pressure  p. 
This  observation  emphasizes  the  need  for  inclusion  of  pressure 
dependent  shear  strength  in  constitutive  modeling  of  poly¬ 
crystalline  TiB2  applicable  to  the  shock  loading  regime.  A  simpler 
model  wherein  yield  strength  Y  and  hence  r  =  ViY  are  constant  for 
shocks  at  or  exceeding  the  HEL,  as  has  been  assumed  elsewhere  as 
an  approximation  for  some  metals  [17],  would  be  highly  inaccurate 
for  the  ceramic  TiB2. 

Predicted  cumulative  inelastic  strain  f  is  shown  versus  shock 
pressure  P  in  Fig.  3.  Below  HI,  no  inelastic  deformation  occurs. 
Between  HI  and  H2  (first  and  second  yield),  the  increase  of  inelastic 
deformation  with  increasing  stress  is  relatively  rapid.  Above  H2,  the 


Fig.  4.  Predicted  (a)  entropy  p  normalized  by  specific  heat 


Fig.  5.  Axial  (Hugoniot)  stress  at  first  yield  point  am  versus  static  lateral  pre-stress  (Tq: 
model  and  experiment  [43]. 

slope  of  f  versus  P  is  less  than  that  below  H2.  Inelastic  deformation 
is  reduced  by  compressive  pre-stress. 

Recovery  and  spall  experiments  [7]  have  demonstrated  that 
micro-cracking,  notably  transgranular  fracture,  dominates  the  in¬ 
elastic  response  at  shock  pressures  between  HI  and  H2,  leading  to  a 
decrease  in  spall  strength.  It  has  also  been  suggested  that 
dislocation-based  slip  may  be  a  primary  inelastic  mechanism  above 
H2  [41].  Shear  stress  at  H2,  on  the  order  of  5  GPa  in  Fig.  2(b),  is 
similar  in  magnitude  to  that  suggested  for  dislocation  glide  resis¬ 
tance  in  alumina  [42].  Recent  shock  experiments  on  titanium 
diboride  [43],  in  which  compressive  pre-stress  was  applied  by 
shrink  fitting  a  metal  ring  around  a  cylindrical  specimen,  have 
demonstrated  an  increase  in  spall  strength  with  compressive  pre¬ 
stress.  This  result  is  consistent  with  model  predictions:  gq  sup¬ 
presses  micro-cracking  associated  with  f,  leading  to  higher  spall 
strength  upon  tensile  release. 

Entropy  and  temperature  predictions  are  given  in  Fig.  4(a)  and 
(b),  respectively.  Temperature  rise  due  to  the  elastic  precursor  is 
approximately  3  K  for  titanium  diboride.  Slopes  of  entropy  and 
temperature  rise  versus  volume  change  increase  with  increasing 
compression  as  contributions  from  inelasticity  and  thermoelastic 
coupling  both  become  large.  Temperature  doubles  to  590  K  at  Vj 
Vo  =  0.915  in  the  material  without  pre-stress.  Pre-stress  results  in 
predicted  increases  in  both  entropy  and  temperature. 

Fig.  5  compares  model  predictions  of  axial  stress  P  at  initial  yield 
(dHi)  with  experimental  data  [43]  at  various  levels  of  pre-stress  gq. 
In  these  experiments  on  titanium  diboride  [43],  compressive  lateral 
pre-stress  was  applied  by  shrink  fitting  a  metallic  ring  around  a 


and  (b)  temperature  d  normalized  by  initial  temperature  6o. 
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Table  3 

Model  predictions  versus  pre-stress  (Tq. 


Variable 

Description 

Value 
((To  =  0) 

Value 

((To  =  0.5  GPa) 

Value 

((To  =  1.0  GPa) 

Po 

pre-shocked  mass 
density  [g/cm^] 

4.510 

4.518 

4.526 

O'Hl 

initial  yield  (HEL)  [GPa] 

6.00 

8.47 

9.58 

O'HZ 

second  yield  point  [GPa] 

14.78 

15.66 

16.45 

®E 

precursor  shock 
velocity  [km/s] 

11.77 

11.95 

11.93 

Ve 

precursor  particle 
velocity  [km/s] 

0.113 

0.157 

0.177 

20  GPa  shock 
velocity  [km/s] 

9.209 

8.859 

8.726 

v~ 

20  GPa  particle 
velocity  [km/s] 

0.450 

0.401 

0.377 

cylindrical  ceramic  specimen,  as  noted  already.  The  HEL  of  the 
material  without  pre-stress  in  these  experiments  is  5.6  GPa,  slightly 
lower  than  the  model  value  of  6  GPa  that  follows  from  Ref.  [5]. 
Regardless,  agreement  between  the  present  model  and  experi¬ 
mental  data  shown  in  Fig.  5  is  considered  close,  with  the  model 
predicting  the  correct  trend  and  magnitude  of  increase  in  HEL  with 
increasing  pre-stress.  Errors  in  the  predicted  HEL  compared  with 
experimental  data  (computed  in  the  same  way  as  in  Table  2) 
are  +7%,  -1%,  and  -2%  for  pre-stress  values  of  0,  0.34  GPa,  and 
0.68  GPa,  respectively. 

Table  3  lists  various  solution  variables  for  different  levels  of  pre¬ 
stress.  Compressive  pre-stress  increases  both  the  first  and  second 
yield  points  and  the  particle  velocity  of  the  elastic  precursor,  the 
latter  computed  from  (3.24).  Pre-stress  has  little  effect  on  predicted 
precursor  wave  speed.  Pre-stress  results  in  a  decrease  in  shock 
velocity  and  particle  velocity  from  (3.23)  for  a  20  GPa  plastic  shock. 
Shock  velocity  without  pre-stress  is  in  reasonable  agreement  with 
experiment  [22]. 

5.  Conclusions 

A  finite  strain  model  has  been  developed  for  brittle  solids  and 
applied  to  describe  shock  compression  of  the  polycrystalline 
ceramic  titanium  diboride.  Analytical  solutions  to  the  planar  shock 
problem  have  been  obtained  for  material  with  possible  lateral  pre¬ 
stress.  Two  parameters  associated  with  pressure  dependent  shear 
strength  in  regimes  below  and  above  second  yield  are  fit  to 
experimental  Hugoniot  data.  The  model  and  solution  then  enable 
prediction  of  mean  pressure,  shear  stress,  inelastic  strain,  entropy, 
temperature,  and  effects  of  pre-stress,  all  without  further  calibra¬ 
tion.  Results  closely  follow  experimental  observations,  including  an 
increase  in  the  HEL  and  decrease  in  inelastic  deformation  with 
increasing  compressive  pre-stress. 
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